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Abstract. As a first step of a program aimed to the detec- 
tion of dark matter (or radial variations of M/L) in early- 
type galaxies, we report deep spectroscopic observations 
of the bulge-dominated edge-on SO galaxy NGC 3115, 
made at ESO, La Silla, using EFOSC at the 3.6 m tele- 
scope and EMMI at NTT. Such observations allow mea- 
surements of the rotational velocity out to 1.8a e (effective 
radii) from the galaxy center, where the surface bright- 
ness is \ib — 24 mag arcsec -2 . The rotation curve quickly 
reaches an asymptotic value, (vf) ~ 260 km s , with only 
marginal indication of systematic decline within the range 
of our observations. The line-of-sight velocity dispersion 
has also been measured; it decreases steeply from a rather 
high central value and flattens out ((a) ~ 100 km s -1 ) 
within our observing range (a ^ 1.3a e ). 

Models built on these data and simple dynamical ar- 
guments show that the M/L of NGC 3115 must thus be 
increasing from M/L = 6 (in solar units) in the inner 
regions (~ la e ) to at least M/L > 10 in the outermost 
regions (~ 2a e ). 

Key words: Galaxies: NGC3115; kinematics and dynam- 
ics of; lenticular; structure of; dark matter 



1. Introduction 

Some evidence for the existence of dark matter has now 
been acquired on all scales, from those of dwarf galaxies 
up to the richest galaxy clusters. At the same time infla- 
tionary theories suggest that dark matter might exist even 
at cosmological scales. However, at galactic scales, that is 
at galactocentric distances comparable with the lengths 
of the effective radii a e , such evidence is still patchy (c/. 
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Kormendy & Knapp 1987). It seems rather firm for dwarf 
and spiral galaxies, whereas no direct evidence has yet 
been discovered for early-types (E and SO's). This hap- 
pens primarily for three reasons. In the first place, early- 
type galaxies have little, if any interstellar matter (Knapp 
1987), which traces the rotation curve and, indirectly, the 
dark matter in spirals; the ionized gas in early-type galax- 
ies is mostly confined to the innermost regions, where dark 
matter is unlikely to be dominant, and all this gas has 
a very uncertain geometry, which severely limits the in- 
terpretation of observations. Cooling flows in ellipticals 
provide some support for dark matter at very large radii 
(Fabian, Arnaud & Thomas 1987), but these estimates 
are fraught with uncertainties (Bertin, Pignatclli & Saglia 
1992). Furthermore they apply only to a small and some- 
what peculiar subset of all ellipticals, the brightest cluster 
members, which all happen to be "boxy" galaxies (Bender 
et al. 1989, Capaccioli, Caon & D'Onofrio 1992). For the 
other class of ellipticals, the "disky" E's, which constitute 
~ 50% of all E's (Bender et al. 1989), nothing is known. 



Bulge-dominated edge-on SO's (hereafter referred to 
as BEDOS) provide a convenient laboratory to investi- 
gate the presence of dark matter inside the optical re- 
gions of early-type galaxies. First of all, bulge dominated 
SO's are now believed to be closely related to the class 
of "disky" ellipticals (Capaccioli 1987, 1990, Capaccioli, 
Caon & Rampazzo 1990) for which even the equivocal 
evidence provided by the X-ray coronas is not available. 
Moreover, the existence of dark matter haloes surround- 
ing SO galaxies seems to be demanded by those few such 
objects possessing polar rings (e.g. Sparke 1989 and Arn- 
aboldi et al. 1992). In this respect, the study of BEDOS by 
direct kincmatical mapping along the major axis is greatly 
simplified for the following reasons: 
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1. they have a high v max /a c ratio (Capaccioli 1979, 
Davics et al. 1983), which makes the analysis simpler (and 
more similar to that of spirals); 

2. the geometry is known since, as i ~ 90°, the appar- 
ent axial ratio equals the intrinsic one; 

3. the bulge is unlikely to be strongly triaxial (for 
NGC 3115 and NGC 3379, see Capaccioli et al. 1991); 

4. direct observations provide the photometric charac- 
teristics and the M/L of the disk (and thus its dynamical 
relevance), which allow us to subtract its effects; 

5. thanks to the "cylindrical rotational field" (i.e. very 
slow variation of the velocity perpendicularly to the galaxy 
major axis; cf. Kormendy & Illingworth 1982), we know 
how to obtain the true rotation velocity of the bulge, with- 
out contamination from the disk; 

6. lastly, substantial photometric and spectroscopic 
literature exists, which permits relevant checks and com- 
plements modern observations. 

A kinematical project involving BEDOS may be simpler 
than considering directly elliptical galaxies, but it is not 
altogether simple. In fact, detection of dark matter re- 
quires spectroscopic observations at large distances from 
the centers of galaxies. As a way of comparison, existing 
models of the mass distribution of the Galaxy (Caldwell 
& Ostriker 1981, and Bahcall, Schmidt & Soneira 1982) 
show that, within the solar radius (8.5 kpc, corresponding 
to ~ 2.5 disk scale-lengths), the mass is roughly equally 
distributed between spheroid, disk, and dark halo. For a 
similar distribution of matter in early-type galaxies, one 
should find ~ 1/2 of the total mass in dark matter within 
~ 1.8 effective radii a e , the radial range corresponding to 
the same fraction of the total light in the Milky Way. Spec- 
troscopic observations of stellar lines at these distances 
from the center of the galaxy are hampered by the faint- 
ness of the continuum, which is notoriously already faint 
enough at 1 a e : fi e ~ 22.5 ±2 B-mag arcsec -2 for E galax- 
ies and SO bulges with total luminosity Mb i5 — 18 (Ca- 
paccioli & Caon 1991; their Fig. 1). 

Kinematical data out to galactocentric distances of 
~ 2a e have now been obtained for the first time as part 
of a program started at ESO, La Silla, in 1989, and aimed 
at investigating the dependence of the dynamical behav- 
ior of galactic disks on the bulge-to-disk ratios, and the 
presence and shape of dark haloes in early-type galax- 
ies (Cappellaro, Capaccioli & Held 1989, 1990) through 
observations of "deep" rotation and velocity dispersion 
curves in BEDOS and polar ring bulges (Arnaboldi et 
al. 1992). This program has been made possible by the 
unique capabilities of the ESO Faint Object Spectrograph 
and Camera (EFOSC; Melnick, Dekker & D'Odorico 1989) 
attached to the Cassegrain focus of the ESO 3.6 m tele- 
scope. The first target galaxies, NGC 2310, NGC 3115, 
and NGC 4179, were chosen among the few BEDOS with 
fairly large angular sizes (distances < 25 Mpc) with the 
purpose of minimizing resolution problems. 



Here we concentrate on NGC 3115, leaving the remain- 
ing objects for a following paper. One reason for this choice 
is that much work on this "standard" SO galaxy has al- 
ready been carried out (see Capaccioli, Held & Nieto 1987, 
for an extensive photometric mapping, and §2 for refer- 
ences on kinematics and dynamics). In particular, we al- 
ready know the mass-to-light ratio of the exponential disk 
of this galaxy (Capaccioli, Vietri & Held 1988) and we 
have also shown that the dominating bulge is unlikely to 
be very triaxial (Capaccioli et al. 1991). Another reason is 
our need to demonstrate, with a good example, the relia- 
bility of the EFOSC data, particularly for what concerns 
the velocity dispersion measurements (see §2), taking ad- 
vantage of the higher resolution but more expensive results 
from spectra taken with EMMI at the ESO NTT. 

The plan of the paper is as follows. Section 2 de- 
tails the spectroscopic observations and their reduction, 
with results doubling the radial extension of the pub- 
lished kinematical measurements in this very well studied 
galaxy. Major-axis rotation and velocity dispersion curves 
of NGC 3115 are analyzed theoretically in Sect. 3. Sec- 
tion 4 contains a summary and our conclusions. An Ap- 
pendix discusses the theoretical method employed for the 
dynamical analysis, and expands on a few technical de- 
tails. The relevant data on NGC 3115, mostly taken from 
Capaccioli et al. (1987), are summarized in Table 1. 

2. Spectroscopic observations 

2.1. The EFOSC spectra 

Five long exposure spectra of NGC 3115 (Table 2) were 
obtained with EFOSC at the ESO 3.6 m telescope, and 
recorded with a high resolution RCA CCD (640 x 1024 
pixels; ESO code #8), in binned mode: pixel size = 30/im, 
or 0"675 on the sky. We chose the ESO Orange 150 grism, 
with dispersion of 3.9 A px" 1 from 5000 to 7000 A. The 
3(6 x 1''5 spectrograph slit was aligned with the major 
axis of NGC 3115 (P.A.= 43°5). Although the nucleus of 
the galaxy was set close to one end of the slit, the other 
end still viewed a region of the object as bright as hb — 
24.5, i.e. only ~ 7 times fainter than the average night- 
sky brightness. Therefore we also obtained long exposure 
spectra of the blank night-sky (by offsetting the telescope 
by ~ 20'), in chronological sequence with the astrophysical 
exposures. 

The raw data were pre-processed by the standard pro- 
cedures in MIDAS (ESO Operating Manual No. 1) for bias 
and dark subtraction, and for flat -fielding based on dome 
and dawn-sky exposures at different count levels. The 
complex distortion pattern was mapped using the com- 
parison spectra to derive a line-by-line wavelength cal- 
ibration. We found that the position of the comparison 
lines drifts over the detector during the night, the am- 
plitude of the displacement being of the order of 1 pixel 
(equivalent to ~ 200 km s _1 ). The consequent frame-to- 
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Table 1 

NGC 3115: catalog data 



Parameter 


Symbol 


Value 






^(^ h 09 m 44 s ^1 


Declination ("1950 0~l 


§ 


—07° 28' 30"7 


Morphological type 




S0"(sp) 


Standard isophotal major diameter 


D 25 


8!63±0!27 


Standard axis ratio 


(b/a) M 


0.51 ±0.03 


Standard position angle 


P.A.25 


O 

43:5 ± 1° 


Corrected total magnitude* 




9.65 ±0.1 


Observed radial velocity [km s _1 ]f 


vo 


663 


Adopted distance [Mpc] 


A 


10 


Scale factor [pc] 


1" 


48.5 


Central rotational vel. gradient [km s~ arcsec - ]f 


I dv(a) Ida ] 

V Ja=0 


> 100 


Asymptotic rotational velocity [km s Jj 


Its \ 

\ v f) 




Central velocity dispersion [km s If 


o~(U) 


ozo 


Velocity dispersion at la e [km s _1 ]f 


cr(a e ) 


105 


Spheroid: 






Effective semi— major axis 


a e 


yo ±0 


Effective semi— minor axis 


i. 

Oe 


OO zt 


Corrected total magnitude 




9.71 ± 0.1 


Mass-to-light ratio at la e [solar=l]f 


M/L 


6 


Mass-to-light ratio at 2a e [solar=l]f 


M/L 


> 10 


Disk: 






Fractional luminosity 


ko 


0.06 ±0.01 


Face-on scale length^ 


h 


25!' 5 


Inclination} 


i 


86° 


Mass-to-light ratio [solar=l]J 


M/L 


7 



* From B T = 9.75 (Capaccioli et al. 1987) with A B = 0.1 (Burstein & Heiles 1982). 
f This paper. 

J From Capaccioli, Vietri & Held (1988). 



frame shift was corrected by matching the position of the 
night-sky emissions. The drift acted on the zero point 
of the absolute velocity scale, adding an uncertainty of 
±50 km s -1 {r.m.s.) to the systemic velocity. More impor- 
tantly, it broadened significantly the long exposure spectra 
of the galaxy, thus introducing a difference with the very 
short exposure spectra of the template stars (which did 
not suffer any significant widening by the drift). Reper- 
cussions of such a difference on velocity dispersion mea- 
surements will be considered below (Sect. 2.4). 

A crucial step in the reduction was the night-sky sub- 
traction from galaxy spectra. We used the blank sky expo- 
sures since no part of the slit in our spectra of NGC 3115 
was free from the galaxy signal. The frames corresponding 
to the same side of the galaxy major axis were individu- 
ally sky subtracted, and then averaged with an algorithm 
correcting for cosmic ray events. 

The kincmatical data were analyzed with a package 
(Bender 1990) based on the cross-correlation (CC), at 
each radial bin, of the continuum-subtracted galaxy spec- 



trum with a template. No differences were found by chang- 
ing the template stars (spectral classes from G8 to Kl). 

2.2. B&C and EMMI spectra 

Given the large uncertainty in the zero point of the veloc- 
ity scale, and the off-centering of the galaxy nucleus at one 
end of each spectrum (which, with the steep central gra- 
dient of the rotational velocity, prevented the sliding fit of 
the velocity curve and the search for the systemic velocity 
by the folding technique — see below), and most impor- 
tantly, the extra-broadening of the galaxy spectra with re- 
spect to the templates, we complemented the deep EFOSC 
material with other spectra of NGC 3115 at higher reso- 
lution. 

In March 1990 we took a relatively 'light' spectrum of 
NGC 3115 with the Boiler & Chivens (B&C) spectrograph 
attached to the Cassegrain focus of the ESO 3.6 m tele- 
scope (Table 2), using grating #10 which, in the second or- 
der, gives a dispersion of 59 A mm -1 over the range 4755- 
5640 A. The 6' long slit was opened to 222 ^im ( = 1"6), 
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Table 2 galactocentric distance up to a maximum of Aa = 10"; 

the r.m.s. deviation about the mean value never exceeds 
NGC 3115: journal of observations 10 km s _1 , even in the most distant bin. 



Instrument 


Nucleus 


Exp. time 


Date 




placed at 


[seconds] 




ESO 3.6 m + EFOSC 


SW 


3600 


1989 Feb. 7 


ESO 3.6 m + EFOSC 


SW 


1536 


1989 Feb. 7 


ESO 3.6 m + EFOSC 


SW 


1800 


1989 Feb. 7 


ESO 3.6 m + EFOSC 


NE 


3600 


1989 Feb. 8 


ESO 3.6 m + EFOSC 


NE 


3000 


1989 Feb. 9 


ESO 3.6 m + B&C 


center 


1200 


1990 Mar. 21 


ESO NTT + EMMI 


NE 


3600 


1991 Dec. 30 


ESO NTT + EMMI 


NE 


3x3600 


1991 Dec. 31 



In all cases the slit has been oriented at P. A. = 43. 5 



aligned with major axis and centered on NGC 3115 nu- 
cleus. The spectrum of the galaxy (and of a few template 
stars) was recorded on the same CCD #8 as above, spa- 
tially binned to 30 /um., with resolution FWHM- 1.75 A 
(= 2 pix). This material was preprocessed and reduced by 
the same technique applied to the EFOSC spectra. The 
systemic velocity, obtained by matching the two sides of 
the rotation curve opposite to the center, is vq — 663 ± 5 
km s _1 , in agreement with Rubin, Peterson & Ford (1980). 
We note explicitly that the coordinate of the velocity cen- 
ter coincides with that of the peak of a, in agreement with 
the symmetry of the velocity dispersion. 

The B&C data were not extended enough to solve all 
problems of the EFOSC spectra. Thus four spectra of one 
hour each (Table 2) were taken in December 1991 with 
EMMI, a spectrograph attached to the Nasmyth focus 
of the ESO 3.5m New Technology Telescope (NTT). We 
used grating #6 which gives a dispersion of 28 A mm -1 
over the range 4900-5400 A. The 6' long slit was opened 
to 1"0, aligned with major axis of NGC 3115. The galaxy 
nucleus was always set at ~ 100" from the NE end of the 
slit. The spectra of the galaxy (and of two template stars) 
were recorded on ESO CCD #18, with pixels of 19 /im, 
with resolution FWHM- 1.4A (= 2.6 pix). 

The EMMI spectra were again analyzed like the 
EFOSC spectra, and then co-added. One important dif- 
ference here is that, due to the relatively longer slit, we 
performed the blank sky subtraction using the signal at 
the SW edge of the co-added spectrum. The systemic ve- 
locity from the co-added spectrum is v Q = 647 ±4 km s _1 
(error not including zero point uncertainty). 

2.3. The rotation curve 

Radial velocity measurements for the combined EMMI 
spectra are plotted in Fig. 1, after folding about the galaxy 
center and the systemic velocity. The data were also av- 
eraged over intervals with amplitude increasing with the 



Fig. 1. Radial velocities v(a) measured along the major axis of 
NGC 3115 and relative to the systemic velocity, plotted against 
the distance from the center of the galaxy (counted positive on 
both sides opposite to the center). The sign of the velocities 
along the SW semiaxis (where v(a) < 0) has been changed. 
Coding for symbols is: EFOSC = filled circles, EMMI = filled 
diamonds. The solid line, reproducing the plain interpolation 
of the two data sets, runs quite flat from a = 50" ~ a e /2 to 
~ 2a e . Representative error bars are also shown. 

The off-centered EFOSC spectra were combined using 
the B&C and EMMI radial velocity curves as a reference, 
the arbitrariness of the slide fitting technique being min- 
imized by the characteristic run of v(a): a steep central 
rise followed by a constant plateau. The resulting mean 
velocities for the EFOSC spectra, plotted in Fig. 1, were 
binned over 3" radial intervals from a = 19"6 and out to 
70"5, and over 10" intervals out to the record distance 
of ~ 2a e . The r.m.s. deviations are larger than the in- 
ternal error estimates from the CC method (both for the 
radial velocities and for the velocity dispersion measure- 
ments). As an upper limit to the uncertainty of the in- 
dividual measurements we computed the standard devi- 
ations of the differences Av(a) = vne(o-) — vsw(a) and 
Aer(a) = o~ne{q) — o~sw{ a )i obtaining, in both 
value of ~ 20 km s _1 . 

Visual inspection (Fig. 1) shows that the mean curve 
produced by the EFOSC spectra — which has about the 
same resolution of the B&C and EMMI curves in the cen- 
tral region (Fig. 2) — runs essentially flat from 20" out 
to the last observed point, with (v) = 262 ± 9 km s _1 
(the quoted uncertainty is the r.m.s. dispersion). Even 
if the EMMI data are fully consistent, within the er- 
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Table 3 

NGC 3115: kinematical data [km s _1 ] 



Fig. 2. Comparison of our radial velocity measurements 
(EFOSC = x, EMMI = +, B&C = o) with literature data: 
Illingworth & Schechter (1982) = filled circles, Rubin et al. 
(1980) = filled squares (NE side) and diamonds (SW side). The 
steeper inner gradient shown by our data finds confirmation 
in the high resolution rotation curve plotted in Kormendy's 
(1987); see also Fig. 6. 



rors, with the EFOSC data, they suggest that v(a) de- 
creases slowly with the galactocentric distance as: v(a) = 
(252 - 0.09 x a) km s"\ for a > 50". The plain interpola- 
tion of the EFOSC and EMMI data, listed in Table 3 and 
drawn in Fig. 1 as a solid line, runs flat at (i>/) = 253 ± 9 
km s _1 from a Si a e /2 out to ~ 2a e . 

A comparison of our results with the other modern 
data available in tabular form (Rubin et al. 1980, Illing- 
worth & Schechter 1982) is also given in Fig. 2; note the 
steep central gradient, the change of slope at a ~ 2", and 
the secondary minimum at a ~ 45" of v(a). It is of interest 
to remember that in the range of these observations, from 
a ~ l/5a e to ~ 1/2 a e along the major axis, the disk is 
brighter than the bulge (Fig. 3). With respect to preced- 
ing studies our result doubles the interval of galactocentric 
distances over which v(a) has been found to run constant 
in NGC 3115. 

2.4- The velocity dispersion curve 

Mean values of the velocity dispersion a from the EMMI 
spectra, after folding about the galaxy center and aver- 
aging with the same scheme as for the radial velocities, 
are shown in Fig. 4. The r.m.s. dispersion on a, which 
is very small at o ^ 1.4a e (Act < 5 kms -1 ), blows up 
at larger galactocentric distances (where it is also based 
on one side of the galaxy only). Therefore wc have dis- 
carded all data at a > 130"; the reason for the difference 



a" 


v(a) 


a (a) 


a" 


v(a) 


a (a) 








302 


34 


260 


118 


1 


103 


260 


37 


256 


113 


2 


135 


232 


40 


250 


122 


3 


157 


220 


43 


255 


107 


4 


162 


212 


46 


250 


110 


5 


172 


209 


49 


243 


109 


6 


174 


204 


52 


253 


118 


7 


177 


198 


55 


256 


122 


8 


186 


190 


58 


257 


104 


9 


205 


172 


61 


249 


107 


10 


211 


165 


64 


250 


107 


11 


217 


170 


67 


252 


104 


12 


224 


167 


70 


266 


97 


13 


236 


160 


80 


251 


101 


14 


244 


153 


90 


252 


104 


15 


247 


153 


100 


238 


104 


16 


255 


144 


110 


261 


116 


17 


249 


147 


120 


250 


110 


18 


252 


151 


130 


238 




19 


254 


144 


140 


249 




20 


269 


132 


150 


250 




22 


266 


129 


160 


257 




25 


265 


124 


170 


274 




28 


262 


119 


180 


257 




31 


256 


123 


190 


250 





Fig. 3. Fractional surface density of light for the disk of 
NGC 3115 along the direction of the major axis. Data have 
been taken from the photometric study of NGC 3115 by Ca- 
paccioli, Held & Nieto (1987). 
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between radial velocity and velocity dispersion curves is 
that the rapidly decreasing signal-to-noise (S/N) ratio af- 
fects more severely a than v measurements. 



corresponding radial velocity curve). The individual de- 
terminations of a (a) were co- added and averaged using 
the same scheme as adopted for v(a). The resulting val- 
ues cr* (a) were still affected by this convolution, induced 
by the drift of the spectra across the detector during long 
exposures. We postulate that the corrected values a (a) 
depend jon the convolved velocity dispersions through the 
relationa 



1/2 



(1) 



where s is a constant to be found by comparison with 
the EMMI and B&C velocity dispersion curves. With 
s = 125 ± 5 kms -1 , corresponding to a 20% increase 
in the instrumental dispersion, we obtained an excellent 
agreement (Fig. 4 and 5) with the EMMI and B&C curves 
as well as with the less resolved data of Illingworth and 
Schechter (1982). Note that the matched values span the 
interval from a ~ 100 to 300 km s" 1 and that even the 
central cusp of a is in agreement with the high resolution 
data plotted in Kormendy (1987) and Kormcndy & Rich- 
stone (1992); see also Fig. 6. Note also that the EMMI 

I velocity dispersion values arc in perfect agreement with 

the EFOSC data over the entire range, proving that the 

Fig. 4. Same as Fig. 1 for the velocity dispersion. large and somewhat empirical correction indeed works. An 

interpolated curve through the two sets of data is sampled 
in Table 3 and plotted in Fig. 4. 



Fig. 5. Same as Fig. 2 for the velocity dispersion. 

The assembly of the measurements from the various 
EFOSC spectra into a final velocity dispersion curve was 
made easy by the sharp central cusp (whose spatial co- 
ordinate was always found to coincide with that of the 



Fig. 6. Comparison of our kinematical measurements 
(filled squares = EFOSC, filled circles = EMMI) with the 
high-resolution observations by Kormendy & Richstone (1992) 
(open circles). The r 1 ^ 4 scale has been chosen to facilitate com- 
parison. 



1 In a preliminary report (Cappellaro et al. 1989, 1990) we 
simply assumed a = a,/ 1.5 at all distances, from comparison 
with Illingworth & Schechter's (1982). 
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In any case, we have analysed the effects of the instru- 
mental resolution Omst and of the signal-to-noise ratio 
for the EFOSC material. Simulations with our reduction 
package show that the systematic error on the measured 
value of the velocity dispersion is ^ 20% for a > 0A<Ji nst , 
i.e. for a > 80 km s _1 in our case. We have also compared 
the CC package with a Fourier quotient (FQ) package de- 
veloped and kindly made available to us by Dr. G. Gal- 
lctta. In agreement with the results of Larsen et al. (1983), 
we find that CC seems superior to FQ in measuring ve- 
locity dispersions, particularly if they are small. The next 
question is the influence of the signal-to-noise ratio; for 
decreasing S/N, the measured cr bj is dominated by the 
noise. Our simulations show that actually, at all a ^ 1.3a e , 
the signal of our spectra is such that we cannot reliably 
measure velocity dispersions which are iS 100 km s _1 . 
Therefore all values of a at a > 1.3a e are conservatively 
discarded. 

A comparison of our measurements with the same lit- 
erature data of Fig. 2 is provided in Fig. 5. Furthermore, 
comparison with Kormendy & Richstone (1992) is pro- 
vided in Fig. 6. 

3. Theoretical models 

The analysis of the data presented above requires fitting 
the light profiles and the rotation velocity curve with an- 
alytical forms. The photometry (Capaccioli et al. 1987) 
was fitted to the projection of a Jaffe law which, for ease 
of computation of the potential gradients, was modified as 
follows: 



p{m 2 ) = p Q 



where 



m 2 {m 2 + Rj) 



m 2 = R 2 + — 



(2) 



(3) 



with c the intrinsic flattening of the bulge, and Rj a pa- 
rameter to be fitted. This form of the density distribution 
has the same asymptotic behavior as Jaffe's law (Jaffe 
1983) for small and large m's, but leads to simpler po- 
tential gradients. Figure 7 shows the fit to the observed 
minor-axis light profile, which represents the pure bulge 
since NGC 3115 is close to edge-on view. The reduction to 
the major axis, taking into account the variable flattening 
of the bulge isophotes (Table 8 in Capaccioli et al. 1987), 
is also given; the excess light in the observed profile is en- 
tirely due to the edge-on disk (Capaccioli, Vietri & Held 
1988). The best fitting Jaffe radius for the major axis is 
Rj = 92", almost coincident with the effective radius of 
the bulge (Table 1). 

An analytical approximation to the observed rotation 
velocity curve, 



R 



Vrot 



R + R c 



(4) 



Fig. 7. Luminosity and flattening profiles and relative growth 
curve for NGC 3115 (from Capaccioli et al. 1987). The open 
circles reproduce the observed light profiles along the major 
and the minor axis; neither of them has been corrected for 
Galactic absorption (Ab = 0.1, after Burstein & Heiles 1982). 
The solid lines are the Jaffe-model fit to the minor axis and 
its reduction to the major axis (Rj = 92"), done using the 
ellipticity profile of the bulge. The excess of the observed major 
axis light profile over the model is due to the presence of the 
exponential disk (Capaccioli, Vietri & Held 1988). The dotted 
line is the radial run of the apparent ellipticity of the bulge; the 
dashed line gives the relative growth of the luminosity, L/Lt- 



is shown in Fig. 8, with = 260 km s _1 , and R c = 2"5. 

The question we are interested in is whether NGC 3115 
has dark matter within the radius for which observational 
data are available. Thus, in a first approximation, we build 
self-gravitating, constant M/L models, which must satisfy 
the constraints provided by the photometric and spectro- 
scopic observations. 

Models of ellipsoidal systems based on Jeans' equations 
(Binncy & Tremaine 1987) can be built provided that ex- 
tra assumptions are introduced to close the system. One 
reasonable such assumption, discussed by Binney, Davies 
& Illingworth (1990; hereafter BDI), is that the distribu- 
tion function depends only on the two classical integrals 
of motion in an axially symmetric potential, the energy 
E and the ^-component of the angular momentum L z . 
Under such circumstances, the velocity dispersion along 
the z-axis equals that along the i?-axis, thus closing the 
system of Jeans' equations. A variation of the method of 
BDI is discussed in the Appendix. It shall be borne in 
mind that all of our models depend also on the accuracy of 
the assumption that the galaxy is oblate. It is well-known 
that NGC 3115 has a nearly perfectly edge-on disk (Ca- 
paccioli et al. 1988); this, and an argument based on its 
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resemblance to NGC 3379 (Capaccioli et al. 1991), lead 
us to expect that such an assumption should be roughly 
correct. In this case, the apparent flattening of NGC 3115 
coincides with its intrinsic flattening, and can immediately 
be inserted into Eq. |^ and in the computation of the po- 
tential gradients. 

The effect of the disk is included: it has been modelled 
as an exponential disk with zero thickness, scale-length 
r = 25". 5 and M/L = 7 (Capaccioli, Vietri & Held 1988). 
Comparison with models without disk shows that its dy- 
namical influence is negligible at all radii, consistent with 
its accounting for only ~ 6% of the total light. For galax- 
ies with massive haloes, the presence of a disk can only 
make matters worse: in fact, the presence of the disk re- 
duces the bulge contribution in the central parts, and thus 
also reduces the predicted rotation and velocity dispersion 
curves in the outer parts, where the disk contribution has 
died away. 

The main result of this analysis is shown in Fig. 9. 
Here, the observed velocity dispersion along the line of 
sight is compared to the prediction of the model. Two dif- 
ferent curves are presented; the dashed one includes a cen- 
tral mass concentration (black hole as in Kormendy 1987, 
where the relevant parameters are given), with which the 
fit to the kinematical data in the innermost regions of 
the galaxy is improved. The presence of this black hole, 
whose dynamical importance in the outermost regions is 
of course negligible, still indirectly affects our modeling, 
since it reduces the amount of matter needed to fit the 
observational data, and thus reduces the M/L of the visi- 
ble matter, assumed independent of radius. In fact we see 
that the model without the black hole fits the outer data 
somewhat better than the model with the black hole. We 
obtain M/L = 6.5 in solar units for the model without 
black hole, and M/L = 6 for the model with black hole. 

However, the most interesting feature of Fig. 9 is that 
none of the models can follow the roughly flat behav- 
ior of the observed velocity dispersion for R 0.5 Rj. 
It is easy to understand that, given the above assump- 
tions, our data must lead to a physically unacceptable 
model. In fact, we have no evidence for a decrease in ei- 
ther the rotation velocity or velocity dispersion, while in 
an isotropic self-gravitating Jaffe model both must de- 
crease over the radial range observed by us. The rotational 
velocity (= (—Rdtft/dR) 1 / 2 ), which of course is an up- 
per limit to the observed one because of the contribution 
of the velocity dispersion to the model's support against 
self-gravity, is shown in Fig. 8 (dashed line) for the as- 
sumed modified Jaffe model with a run of ellipticity as in 
Table 8 of Capaccioli et al. (1987). It can be easily seen 
that such rotational velocity decreases, contrary to our 
evidence. Thus our data alone are sufficient to exclude a 
self-gravitating model with constant M/L. 

On the basis of our models it is less easy to estimate 
the magnitude of the mass discrepancy in the outer parts. 
Fig. 9 shows that both curves predict Oy — for observ- 



Fig. 8. Rotation curve for NGC 3115. The data points are 
those of Fig. 1 and 2. The solid line plots the approximation 
of Eq. 4, with R c — 2"5. The dashed line is the projected 
rotation velocity for a Jaffe model with zero velocity dispersion 
and with the run of ellipticity variations taken from Fig. 13 of 
Capaccioli et al. (1987). 



Fig. 9. Velocity dispersion along the light of sight measured 
along the major axis Data are from Fig. 4 and 5. The solid 
line is a fit with M/L = 6.5 in solar units, and no central 
black hole; the dashed line is another fit with M/L = 6 and 
a central black hole with parameters taken from Kormendy 
(1987). Velocity units are km s ; radii are in units of the Jaffe 
radius Rj = 92". 
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able radii, R ~ O.SRj for both models. Clearly, such re- 
sults depend on our assumptions, in particular on hav- 
ing assumed the galaxy to be oblate and a 2 = a\. ft 
may be possible, by relaxing such constraints, to build a 
better-fitting model. Still we would like to point out that, 
from our data, there is a minimum M/L for any model 
of this galaxy. In fact, we know that, for R = 180", it 
is v rot w 260 km s _1 . A minimum M/L model will have 
from there on o~ z — = uji = ay = 0; then the minimum 
M/L can be estimated by equating, at that galactocen- 
tric distance, centrifugal and gravitational accelerations. 
We obtain (M/L) m i n = 10 in solar units. It is well known 
that the inner rotation curve (R < 0.5Rj) can be fitted 
with M/L = 6 (Illingworth & Schechter 1982, Rubin et al. 
1980), so that we now find a discrepancy of at least a factor 
of 2 in M/L while moving from R ~ O.bRj to R ~ 2Rj. It 
is easy to understand this from Fig. 7, which shows that 
L changes by only 20% from R = 90" to R = 180", while 
the minimum mass M m i n (R) oc Rv 2 ot doubles over the 
same range, for our flat rotation curve (Fig. 8). 



3. There is a simple argument to show that the disk is 
not only unimportant from the photometric point of view, 
but also from the dynamical point of view. 
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Appendix 

We summarize here the method used to analyze the kine- 
matical data, which is similar to that developed by BDI. 
We model the galaxy as a modified Jaffe model (MJM, 
from now on), with density constant on self-similar ellip- 
soids 



2 2,2, 



(5) 



4. Summary and conclusions 

This is the first paper of a series based on deep major- 
axis spectra of bulge-dominated edge-on SO's, obtained 
by using EFOSC in combination with the ESO 3.6 m tele- 
scope and, in this one case only, EMMI at NTT as a check 
for our unconventional use of EFOSC. Aim of this series 
is to report on direct kinematical evidence of dark mat- 
ter haloes within the optical images of galaxies related to 
"disky" E's, i.e. that half of all normal ellipticals which 
do not offer any other way for probing the occurrence of 
dark haloes. We summarize briefly our conclusions: 

1. We have measured the rotation curve of NGC 3115 
out to a galactocentric distance ~ 2a e , where the surface 
brightness is <~ 20% of that of the night sky. It increases 
steeply in the galaxy core, then remains constant (within 
10%) at (vf) = 260 kms" 1 out to 2a e . Similarly, the 
observed velocity dispersion peaks sharply at the center 
(with a maximum of 325 km s _1 ), decreasing outwards till 
it flattens out at the level of <~ 100 km s" 1 ; reliable mea- 
surements are limited to a < 1.3a e . There is a marginal 
indication from the EMMI spectra that a decreases for 
a > 2'. 

2. Our simple minded models, assuming oblateness, 
constant M/L, and a 2 = crfj, predict vanishing velocity 
dispersions along the line of sight, well within our observed 
range. We have argued in a qualitative way that fitting our 
data with M/L=const, is unlikely, even when the afore- 
mentioned assumptions are relaxed. Also, we have shown 
that a simple argument implies that, within our observed 
range, M/L increases by at least a factor of <~ 2, from 
M/L = 6 in the center to M/L = 10 at our outermost 
radius (1.8a e ). 



(c is the true axial ratio) and density profile 



2tt 2 c (i? 2 + Ri)(m 2 + R 2 )(m 2 + R 2 ) 



(6) 



where the total mass M = 1 and i?2 2> Ri are two con- 
stants. The introduction of a second scale parameter (i?i , 
much smaller than any observed radius) is used to simplify 
the computation of the gravitational potential. 

The model that we are trying to fit to the observa- 
tional data must satisfy Jeans' equations. Now we assume 
that the system that emits the light is oblate and self- 
gravitating (i.e. it has constant mass-to-light ratio), and 
has velocity dispersion v 2 along the z-axis equal to that 
along the i?-axis, v 2 R . The first two assumptions are obvi- 
ous, the third one is discussed by BDI. Then Jeans' equa- 
tions for this model in cylindrical coordinates are 



d(p° 2 R ) 
OR 



d(pa 2 z ) +p d(l 



'4 , d($ + (t> ext ) 



R 



dR 



= 



'Jext) 



dz 



dz 



= 



(7) 



(8) 



Here <p is the gravitational potential sclf-consistcntly en- 
gendered by the light -emitting system, and is given for 
reference below in explicit form in the case in which c 
is independent of radius, and <j>ext 

is the disk potential, 
computed for an infinitely thin exponential disk (Binncy 
& Tremaine 1987), with the parameters described in Sec- 
tion 4. Clearly, since p and <f> are supposed known, these 
two coupled equations become a system of two unknowns, 
o~ 2 = a\ and u?, which can be readily solved numerically. 
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The total velocity dispersion thus determined in- 
cludes contributions from both the rotational velocity and 
the true velocity dispersion: 



v l = v rot + a \ 



(9) 



On the other hand, observations give us v 2 ot , which can be 
subtracted from v 2 ^ to yield a 2 ^. We now have all the prin- 
cipal components of the velocity dispersion tensor, and 
can thus project them (see for instance Binney & Mamon 
1982) to give a v , the velocity dispersion along the galaxy 
major axis, and this can be compared with the observed 
quantity <r 2 bs . However, we still have the extra freedom of 
an unknown multiplicative constant, essentially the M/L 
of the galaxy. In keeping with the spirit of the model, 
which is self-gravitating and thus ought to reproduce the 
kinematics of the inner regions of the galaxy, we fix the 
M/L by demanding that the theoretically determined o~y 
fits a 2 bs in the inner regions of the galaxy. 

Now our model has no free parameters left, and ought 
to reproduce a 2 bs in the whole observable region. 

The main difference between our method and that of 
BDI lies in the fact that we treat the surface photome- 
try as if the light profile where of the Jaffe type every- 
where, while they use a Lucy algorithm to deproject the 
two-dimensional image (or close to this: they have light- 
emission profiles along a number of PAs). We compensate 
for this by taking for the axial ratio c not a constant, but 
the actual observed run with R. Since the determination 
of c in standard packages essentially requires observations 
along 3 PA's (we discard the term oc cos 49), this com- 
pares reasonably well with BDI, who had data along 7 
PA's. Also, since NGC 3115 is seen edge-on, the axial ra- 
tio c is not a free parameter for us, since c = c f, s . Lastly, 
BDI do not subtract the observed velocity rotation curve 
like we do. But since in this case the data for a 2 bs and 
v rot are coextensive, the method adopted here seemed the 
most appropriate one. 

We give here for reference the gradients of the gravi- 
tational potential engendered by the flattened MJM, for 
the case of constant axial ratio c. In these formulae, some 
simplification is possible through the fact that the density 
(Eq. ^) can be rewritten as 



1 



27r 2 c (Ra+Ri) 



1 



1 



(10) 



Such potential is available through well-known implicit 
formulae {e.g. Binney & Tremaine 1987). We define R 2 = 
x 2 +y 2 , take G = 1, and omit the algebra, to find 



dj>_ _ 
8R ~ 

where 



2R 



7T {R 2 - Rl) 



X(R 1 )-X(R 2 ) 



(11) 



X(Ri) = -^ «?(T 1 ,T 2 ,l)+ ff (T 2 ,l,T 1 )+ ff (l,T 1 ,T 2 ) (12) 

Ri 



Here we have defined n < t 2 as the moduli of the solutions 
of 



R 2 t 2 + (l + c 2 )R 2 + z 2 +R 2 t + R 2 c 2 + z 2 + R 2 c 2 = 0(13) 

(which, of course, depend on i?j), and 
(1-c 2 ) 1 / 2 



g(Tl,T 2 ,X = 

(T 2 - Tl)(X - Tl) 

Analogously, we find 
86 2z 



\ - arctan ( Tl - C c2 )i/ 2 



(14) 



dz n(R 2 - i?i) 
with 

1 



Y(Ri) 



h(r) ee 



Rt(n-T 2 ) 



— arctan ■ 



r(i2i)-r(JZ2) 



h(r 2 ) - h(n) 



(r-O 



2TT72 



(15) 



(16) 



(17) 



( r _c3)l/2 

and the r's are defined as above. 
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